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Abstract 

We propose a set of conventional Bethe Ansatz equations and a corresponding 
expression for the eigenvalues of the transfer matrix for the open spin-^ XXZ quantum 
spin chain with nondiagonal boundary terms, provided that the boundary parameters 
obey a certain linear relation. 



1 Introduction 



Consider the open spin-| XXZ quantum spin chain with nondiagonal boundary terms, de- 
fined by the Hamiltonian P] 

^ N-l 

^ = 2! ^ («+i + ^K+i + coshr/ 

n=l 

r 2k- 

+ sinh 71 coth f _crf H ; f cosh 9_a^ + i sinh ^-crF) 

L ' sinh ^_ ^ ^ 

- coth ^+(7^ - ^.^^ ^ ( cosh g+q-^ + i sinh 6'+cr^) |, (1.1) 

where , o"^ , are the usual Pauh matrices, t] is the bulk anisotropy parameter, , k± , ^± 
are arbitrary boundary parameters, and is the number of spins. This is the prototypical 
integrable quantum spin chain with boundary. It is related to many other models, including 
the sine-Gordon field theory Moreover, this model has applications in various branches 
of physics, including condensed matter and statistical mechanics. 

This model has resisted solution for many years (see, e.g., jH])- The main difficulty is 
that, in contrast to the special case of diagonal boundary terms (i.e., n± = 0) considered in 
jUE], a simple pseudovacuum (reference) state does not exist. For example, the state with 
all spins up (g) is not an eigenstate of the Hamiltonian. 

We recently formulated jH] a method of deriving the Bethe Ansatz solution of integrable 
spin chain (vertex-type) models which does not rely on the existence of a pseudovacuum 
state. In particular, we used this method to solve the model (jl.lj) for the special case 

= e+ = e-, ^+ = ^- = 0, iV = odd. (1.2) 

Here we propose the solution for a more general case. Indeed, in terms of the boundary 
parameters a=p , introduced below in Eq. ()H.24j) . we find an expression for the eigenvalues 
of the transfer matrix corresponding to the Hamiltonian 

^=2(5^ («+i + + cosh ri 

n=l 

-|- sinh 1] ^ coth a;_ tanh P-a^ + csch a_ sech /3_ ( cosh O^a^ + i sinh 9-crf) 

— coth a+ tanh /3+(T^ + csch a+ sech /3+ ( cosh 6'+cr^ + z sinh 6'+cr^) |, (1.3) 
where the boundary parameters are subject to the linear relation 

a- + (3- + a+ + P+ = ±{e--e+)+r]k, (1.4) 
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where k is an even integer if is odd, and is an odd integer if is even. In the recent paper 
jZj, similar results have been obtained by a different approach. 

The outline of this article is as follows. In Section |21 we briefly review the construction of 
the model's transfer matrix, and list some of its important properties. In Sectional we find 
the eigenvalues of the transfer matrix by the three-step procedure formulated in The 
first two steps, which lead to a functional relation for the transfer matrix, are the same as in 
inilHI, except for the introduction of the parameters 9^. The principal new results appear at 
the third step, where we succeed to recast the functional relation in terms of a determinant 
for the more general case ()1.4|) . We conclude with a brief discussion of our results in Section 

m 



2 The transfer matrix 



The fundamental transfer matrix t{u) corresponding to the model is given by [S] 

t{u) = tio K+{u) To(n) K,{u) fo{u) , (2.1) 
where the monodromy matrices are given by 



To(m) = Ron{u) ■ ■ ■ Roi{u) , To(m) = Rio{u) ■ ■ ■ Rm{u) 



(2.2) 



and the R matrix is the solution of the Yang-Baxter equation given by 

\ 

n sinVi n «inVi n D 

R{u) 



( sinh(u + ?7) 

sinhti sinh?7 

sinhr; sinhw 

V sinh(M + r/) / 



(2.3) 



Moreover, the K matrix is the solution of the boundary Yang-Baxter equation jH] given by 



K-{u] 



sinh(^_ + u) K_e^" sinh2'u 
K-e~^~ sinh2M sinh(^_ — u) 



(2.4) 



which evidently depends on three boundary parameters C,- , n~ ,6-- It is related to the 
symmetric matrix (u) used in jHl IE] by a gauge transformation, 

6— =0 



K-(u) = M K-(u) 



(2.5) 
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with 



The matrix K~^{u) is equal to K^{—u — 77) with (^_,k_,6'_) replaced by (^+,k+,6'+). 
Finally, trg denotes trace over the (two-dimensional) "auxiliary space" 0. Further details 
about the construction of this transfer matrix can be found in |H] . 

The transfer matrix constitutes a one-parameter commutative family of matrices 

[t{u),t{v)] = Q. (2.7) 
The Hamiltonian (jl.lj) is related to the first derivative of the transfer matrix, 

+ C2I, (2.8) 



u=0 



where 



1 sinh^ 1] + N cosh^ i] 

4 sinh ^_ sinh ^+ sinh^^"^ r] cosh r] ' 2 cosh 7] 



and I is the identity matrix. The two relations ()2.7j) . ()2.8j) signal that the model is integrable. 
Moreover, it is evident that in order to determine the energy eigenvalues, it suffices to 
determine the eigenvalues of the transfer matrix. 

The transfer matrix has the periodicity property 

t{u + i7T) = t{u) , (2.10) 

as well as crossing symmetry 

t{-u - 7]) = t{u) , (2.11) 
and the asymptotic behavior (for k± 7^ 0) 

^u{2N+4)+r,{N+2) 

t{u) fi;_/t+cosh(6'_ - 6*+) ^j^^ for u ^ 00 . (2.12) 



3 Bethe Ansatz solution 

We now proceed to find an expression for the transfer matrix eigenvalues using the method 
formulated in jS]. This method consists of three main steps: 
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3.1 Step 1: fusion hierarchy 



The first step is to obtain the model's so-called fusion hierarchy. The transfer matrix (j^.lj) 
is actually the first (j = |) member of an infinite hierarchy of commuting transfer matrices 
t^^\u) corresponding to spin-j (i.e., (2j + l)-dimensional) auxihary spaces, j = | , 1 , | , . . .. 
Using the fusion procedure for R jTUl E] and K ^21 UH] matrices, one finds that these 
higher-level transfer matrices obey the relations 



t^'\u) = C2j-i(2m + (2j - 1)7]) [t^'--^\u) t^-2\u + (2j - 1)7]) 

A{u + (2j - 2)7]) C2,_2(2^ + (2j - 2)7]) 

C(2« + 2(2j - 1)7]) ^ ! 



with = I, and j = 1 , | , . . .. The quantum determinant A('u) is given by 

A(m) = — [sinh('u + ?7 + ^_) sinh(M + ?7 — + sinh^(2'u + 27])] 
X [sinh(M + T] + ^+) sinh(u + r] — + sinh^(2M + 27])~\ 
X sinh 2u sinh(2u + Arj) ({u + 7])'^'^ , 

and 



(3.i: 



(3.2) 



OH 
C(«) 



l[au+k7]) 



Coin) 



k=l 



— sinh(M + 7]) sinh(u — 7]) . 



(3.3) 
(3.4) 



These relations are the same as those for the case of symmetric K matrices {6zp = 0) 
We remark that the spin-j matrix K^^ 2j)(^) related to the corresponding matrix with 
= by a generalization of the gauge transformation ()2.5p . 



^ri...2,>(«) 



Ml--- M2j ^(i...2,)(«) ---Mi^ ■ 



9_=0 



(3.5) 



3.2 Step 2: truncation at roots of unity 

The second step is to observe [8J that for anisotropy values 

in 

V = —T^ P= 1,2,... , 



p+1 



(3.6) 



(and hence g = is a root of unity, satisfying g^"*"^ = — 1), the level- transfer matrix can 
be expressed in terms of a transfer matrix of one level lower. 



t(^)(M) = a{u) y'^\u + 7])+ P{u)I 



(3.7) 



4 



The quantities a{u) and P{u) are given by the corresponding expressions (4.31) in i^, except 
with cr=p(M) — > e''P^^'*^'Fcr=p(M) and p^{u) e^^^^^^^^^ p^{u), as a consequence of ()3.5|1 . 

This resuh provides an example of McCoy's dictum "Comphcated is simple" ^3] . Indeed, 
the essential point of this step is to exploit the higher symmetry which occurs at roots of 
unity to help solve the model. 

Combining the fusion hierarchy ()3.H) and the truncation identity (or "closing relation") 
fl3.7|) for the 7] values ()3.6|) . we arrive at a functional relation for the fundamental transfer 
matrix t(u) = t^^\u) (and hence, for the corresponding eigenvalues A(m)) of order p + 1 

m- 

A{u)A{u + //)... A{u + p7]) 

- 5{u - r])A{u + r])A{u + 2r]) . . . A{u + {p - 1)?]) 

- S{u)A{u + 2T])A{u + 3r]) . . .A{u + pr]) 

- 5{u + r])A{u)A{u + 3r])A{u + 4:1])... A{u + prj) 

- 5{u + 2r])A{u)A{u + ri)A{u + 4?]) . . . A{u + pr]) - . . . 

- 5{u + (p- l)r])A{u)A{u + //)... A{u + {p - 2)r]) 

+ ... = /(«), (3.8) 

where S{u) is defined by 

C{2u + 2r]) 

Moreover, the function f{u) is given by 

(-1)^(^+1) . cosh^((p+l)n + f6) 

= ^T^^iTvTT^ smh {{p+iy 



22p{N+i) ^ ^ cosh'((p+ l)n) 

X |n(M ; ^_ , K_) n{u ; , k+) + n{u ; , k^) n{u ; ^+ , k+) 

+2(-l)^(-K_K+)P+isinh^(2(p+ l)M)cosh((p+ 1)(^_ - 0+))| , (3.10) 

where e = 2frac(p/2) equals if p is even, and equals 1 if p is odd; and the function n{u ; ^ , k) 
is defined by 

m 

e,/€) = sinh((p+l)(e + M))+ J2 Cp,/ K''sinh((p+l)u+(p + l-2Z)0 , (3.11) 



1=1 
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with 



Cp,l 



{P 



k=0 



For instance, for the case p = 3, the functional relation is given by 
A(m)A(m + r])A{u + 2r])A{u + 3r]) - S{u - r])A{u + r])A{u + 2r]) 

-6{u)A{u + 2r])A{u + 3r]) - 6{u + r])A{u)A{u + 3r]) - 5{u + 2r])A{u)A{u + t]) 

+6{u)6{u + 27]) + 6{u - r/)5(n + r]) = f{u) . (3.12) 



3.3 Step 3: determinant representation 

Following the strategy used in [T3] to solve RSOS models, the third and final step is to 
rewrite the functional relation as the determinant of a {p + 1) x (p + 1) matrix. Let us 
assume that this matrix has the same form as the one for the diagonal case {k± = 0) and 
for the case (|1.2j) . That is, we assume the functional relation can be cast in the form jH] 



det 



Ao 







... 








-ho 


-h 


Ai - 


-K 


... 














-h2 


A2 - 


... 




















... 




Ap-1 


-K- 










... 





-hp 


Ap 




hk = h{u 


+ -qk) 


K = h' 


(m + ijk), 










h\u) 


= h{—u 


-2r/), 







(3.13) 



(3.14) 

and the function h{u) is yet to be determined. We find that the functional relation ()3.8|1 can 
indeed be recast in the form ()3.13p . provided that h{u) satisfies the three conditions 



h{u + m) 
h{u + rj)h{—u — 7]) 



h{u) 
6{u) 



(3.15) 
(3.16) 

(3.17) 



JJ h{u + jri) + JJ h{-u - jr]) 

j=0 j=0 

The results |6j for h{u) in the diagonal case and in the case ()1.2|) suggest that, in general, 
h{u) has form 



, / X • 1 2Ar/ .smh(2u + 2r]) . . , , 

h{u) = - smh (m + m^-TTT, — \ — r9Au)g^{u) 

smh(2u + ri) 



(3.18) 
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where the functions g^{u) contain all the dependence on the boundary parameters. 
Then second condition ()3.16|) together with ()3.2|) . ()3.9p and ()3.18|) imply that 

g^{u)g^{u)g^{—u)g+{—u) = (sinh^ m — sinh^ ^_ + sinh^ 2^) 

X (sinh^ u - sinh^ ^+ + sinh^ 2u) . (3.19) 
This suggests that g^{u) obey the functional equation 

gzf{u)gzf:{—u) = — (sinh^ u — sinh^ + sinh^ 2u) . (3.20) 
Assuming that the functions g^{u) are given by 

g^{u) = 2nzs^ sinh(M + a^) cosh(u + [3^) , (3-21) 
then ()3.20|) implies that the parameters Oip, /^^p obey 

sinh^ cosh^ 13^3- = —5- sinh^ C,^ , cosh^ sinh^ /5=r = — 5- cosh^ ^zp . (3.22) 

A similar reparametrization appears in [3 [7] . Below we shall argue that ()3.21|) is essentially 
the unique solution of (|3.19p . 

The third condition (j3.17j) together with (j3.1(J|) and (j3.18j) imply that ^ 

P V 

Wg-{u + jri)g+{u + JT]) + \Yg^{-u - jri)g+{-u - jri) 

j=0 j=0 

n{u ; ^_ , K_) n{u ; , k+) + n{u ; , k_) ^(m ; , k+) 



22p 

+2(-l)^(-fi:_K+)P+^ sinh2(2(p + 1)m) cosh((p + 1)(0_ - ^+))| , (3.23) 



where n{u ; ^ , k) is given by Eq. (j3.11|) . We find that this requirement can be satisfied for 
p = odd, with gif{u) given by ()3.21|) and ^ 

sinh a_ cosh /3_ = sinh f _ , cosh a_ sinh /3_ = cosh , 

^ 2k_ ^ ' ^ 2k_ ^ ' 

sinh a+ cosh /5+ = sinh ^+ , cosh a+ sinh /?+ = cosh ^+ , (3.24) 

2k+ 2k+ 

provided that the various parameters obey the linear constraint 

a_ + P_ + a+ + P+ = ±{9_-9+)+r]k, (3.25) 

^Note that the right-hand-side of (|3.23l) depends on N only through its parity (—1)^. 
^The requirement of including the special case (|1.2|l . which corresponds to a- ~ —ci+, P- = fc = 0, 

helps to resolve the sign ambiguity in passing from H3.22|l to 13.24|l . 
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where k is an even integer if is odd, and is an odd integer if is even. 

In short, the functional relations can be cast in the determinant form ()3.13p for p = odd 
with 

... . ,2Nf , s sinh(2M + 2r/) 

h{u) = -smh [u + r]) . 

smh(2M + r]) 

X 4k_/«_|_ sinh(M + q;_) cosh(u + sinh(u + a+) cosh(M + , (3.26) 

where a=p, /5=p are defined by ()3.24|) and satisfy the constraint ()3.25|) . 

We now proceed as in OEl, and assume that the matrix in ()3.13jl has the null vector 
iQo,Qi,... ,Qp). That is, 

AoQo - h'_^Qi - hoQp = , 
-hQk^i + AkQk - K^iQk+i = 0, k = l,...,p-l, 

-h'p_^Qo - hpQp^i + ApQp = 0. (3.27) 

We make the Ansatz Qk = Q{u + rjk), where Q{u) is given by 

M 

Q{u) = Y\ sinh(M — Uj) sinh(M + uj + i]) , (3.28) 
j=i 

which has the crossing symmetry Q{u) = Q{—u — r]). The zeros Uj of Q{u) are still to be 
determined. Eqs. fj3.27|) and ()3.14|) imply that the eigenvalues are given by 

AM = /.W%-^ + ft(-«-r,)%i^?l. (3.29) 

Q[u) Q[u) 

We verify that this result is consistent with both the periodicity ()2.10|) and crossing 1)2.111) 
properties of the transfer matrix. The requirement that A(m) be analytic at m = Uj yields 
the Bethe Ansatz equations 

Ku,) _ Q{u,+rj) ^ = i__^M. (3.30) 



h{-Uj - T]) Q{uj -1])' 

The asymptotic behavior ()2.12p . together with the result ()3.29|) for the eigenvalues and the 
constraint ()3.25j) . imply that the number M of Bethe roots is given by 

M = ^{N-l + k), (3.31) 

where k is the integer appearing in ()3.25p . We leave to a future investigation the interesting 
question of determining further restrictions on the value of k, which presumably is related 
to the question of completeness. 
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We now argue that ()3.21|) is essentially the unique solution of fj3.19|) . Indeed, if gz^^u) are 
also solutions of ()3.19j) . then 



9~{u)g+{u) = g^{u)g+{u)(t){u) , (3.32) 
where gzf{u) are given by ()3.2H) . and satisfies 

</)(m)(/)(-m) = 1 . (3.33) 
The periodicity condition ()3.15jl implies that (f){u) has the same periodicity 

0(m + ivr) = 0(m) . (3.34) 

We infer from (HO!?!) and (HOI) that (f){u) is a CDD-like factor 

-p-r sinh{u + Vj) 

0(w)=||— TT -T- (3-35) 

smh(M — Vj) 

The requirement that A{u) be analytic then restricts (f){u) to the form 
q{u - Vj) 



[u] 



q{u) 



where q{u) = Y\ sinh(M — Vj) sinh(M + r] + Vj) . (3.36) 



This is equivalent to having additional Bethe roots, which can be included in Q{u) ()3.28j) . 

Although the above results for the eigenvalues ()3.2(ij) . ()3.29|1 . ()3.3()j) have been obtained 
under the assumption that r] is restricted to the values ()3.6|) with p odd, we expect that these 
results remain valid for generic values of t]. Indeed, we have explicitly verified that these 
expressions reproduce the correct eigenvalues for = (with M = 0) and = 1 (with 
M = 1) for arbitrary 77. 



4 Conclusion 

Our proposed expression for the eigenvalues A(m) of the transfer matrix ()2.1|) corresponding 
to the Hamiltonian ()1.3|1 is given by ()3.29|1 . where h{u) and Q{u) are given by ()3.2(ij) and 
(I3.28|l . respectively; the Bethe Ansatz equations are given by ()3.3()|1 . with M given by ()3.31|1 : 
and the parameters a:^, P^p (which are related to , by ()3.24|) ) must satisfy the constraint 

It remains an open question whether a solution with Bethe Ansatz equations of the 
"conventional" form ()3.30j) can be found which does not require a constraint among the 
boundary parameters. (Although the solution proposed in [Hj does not require any constraint 
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among the boundary parameters, it holds only for the rj values ()3.6p . and the Bethe Ansatz 
equations are not of the conventional form.) 

For the special case (jl.2j) . an analysis of the thermodynamic {N oo) limit and an 
extension to higher- dimensional representations has recently been given in jTHj. For the 
more general case discussed here, it should now be possible to address such questions, and 
also to find generalizations to higher rank-algebras, both for the trigonometric and elliptic 
cases. 
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